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Abstract 

We consider a mesoscopic mechanism of the exchange interaction in a system 
of alternating ferromagnetic/nonmagnetic metaUic layers. In the case of small 
mesoscopic samples the sign and the amplitude of the the exchange interac- 
tion energy turn out to be random sample-specific quantities. They can be 
changed by applying an external magnetic field, by attaching to the system 
superconducting electrodes with different phases of the superconducting order 
parameter and by changing the chemical potential of electrons in the metal 
with the help of a gate. In the case of square or cubic geometries of the non- 
magnetic layer at low temperature the variance of the exchange energy turns 

out to be sample size independent. 
Suggested PACS index category: 05.20-y, 82.20-w 
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The investigation of oscillations of the exchange energy as a function of the nonmagnetic 
layers' thickness and the "giant magnetoresistance" in the system of alternating ferromag- 
netic/nonmagnetic metallic layers has both scientific and technological interest 

The origin of the exchange interaction between the ferromagnetic (F) layers is the 
Ruderman- Kittel type interaction: The interaction of itinerant nonmagnetic metal elec- 
trons with localized "f" or "d" electrons in the ferromagnets induces a spin polarization 
in the nonmagnetic metallic (N) layers. This magnetization, in turn, creates the effective 
interaction between two localized spins in different ferromagnetic layers with the energy 

E{R,,R2)^Iij{Ri,R2)SlS]. (1) 

12 ~* 

Here S'j ' are components of localized spins in the F layers and Ri^2 are their coordinates. 
We employ the simplest model where conduction "s" -electrons interact with locahzed "f or 
"d" electrons in the F layers via a contact interaction with energy 

AY,S{n-R)skS. (2) 

k 

Here and Sk are coordinates and spins of conduction electrons in the metal which are 
labeled by the index k and A is the interaction constant. In the case where two magnetic 
ions are embedded in a 3-d nonmagnetic metal Eqs.1,2 lead to a well known form of the 
exchange energy [^^1 

/.,(4,^2) = /o^^^^^^i^^^. (3) 

^ {PF\R,-R2\r 

Here Iq = || is the interaction energy of adjacent localized spins in ferromagnets, which 
is of order of ferromagnet's critical temperature; Ep, pp are the Fermi energy and the Fermi 
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momentum, respectively, and n = ^ is the concentration of electrons. Following to [^'''''^'^1 
we will use the approximation where the total exchange interaction energy E between the 
magnetic moments in the F layers is a sum of E{Ri,R2) over coordinates ^i, R2 of the 
localized spins in the ferromagnetic layers. 

Y.^ Iij(Ri, 4)5f = Iijnun2j (4) 

Ri,R2 



Here riu and n2j are components of unit vectors Hi, n2 parallel to magnetizations in the first 
and the second layers respectively. We assume that the exchange energy inside ferromagnets 
is large enough and fluctuations of spin magnetization direction along the F layers can be 
neglected. As a result, in the case of clean N layers (L <^ Z), the value of lij and the relative 
orientation of magnetizations of different F layers are oscillating functions of L. Here 
L and / are the thickness and the elastic electron mean free path the in N layer, respec- 
tively. The same effects takes place in the system of alternating ferromagnetic/degenerated 
semiconductor layers. 

Until now, however, both experimental and theoretical studies of this phenomenon have 
been restricted to the consideration of the infinite dimensions of both F and N layers and 
clean N layers. In this article we discuss the opposite case of small sample sizes and disor- 
dered N layers where the mesoscopic effects determine both the exchange interaction between 
the ferromagnetic layers and the conductance of the system. Let us consider, for example, 
a system of two F layers of sizes ^1,^2, -^3 divided by a disordered N layer with L > I (See 
Fig.l). In this case the amplitude of the Ruderman-Kittel exchange interaction between a 
couple of spins located inside different F-layers {Iij{Ri, R2)) ~ exp(— l-^iy-^^l ) exp(— j), 
averaged over realizations of a random potential (and consequently the total average ex- 
change energy (E) between F layers) is exponentially small and can be neglected. I^^^ Here 
brackets () correspond to averaging over realizations of a random scattering potential. On 
the other hand, it is well known [i^-ie] ^^^^ exponential decay of the average {Iij{Ri, R2)) 
is the consequence of randomization of the sign of Iij{Ri, R2) and that the typical amplitude 
of the interaction y^((ii~(^i7^2))^ ~ |^i — ^2!"*^ decreases with distance in the same way 
as in the pure case. Here d is the dimensionahty of space. Therefore, in the case L » I 
the amplitude of the Ruderman-Kittel interaction between the ferromagnetic layers must be 
controlled by the mesoscopic effects. 

In this article we discuss the typical magnitude of the " mesoscopic part" of the exchange 
energy between the ferromagnetic layers. We show that it is determined by long range 
correlations between Iij{Ri,R2) and Ik,i{Rz^Ri) which survive even on distances |^i — 
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RslA-^i — Ra] ^ I- As a result, for example, in the case of square or cubic geometries of 
the N layer and at low temperatures the variance of the exchange energy turns out to be 
sample size independent. 

Let us start with the case where the length L2 of the F-layers is relatively short and 
one can neglect the fluctuations of the orientation of magnetizations along the F layers. 
To find the relative angle 0{ni,n2) between magnetization's angles in different layers one 
has to calculate the sign and the amplitude of the quantity /y. In the case where metallic 
region is disordered, lij is a random sample-specific quantity, which can be characterized 
by its average and moments. One can calculate the correlation function (SlijSlki) with 
the help of the standard diagram technique I^^l for averaging over realizations of random 
potential. Here 5Iij = lij — {lij). The diagrams which contribute to {5Iij6Iki) to lowest 
order in parameter ^ <C 1 are shown in Fig. 2, where solid lines correspond to electron 
Green functions in Matsubara representation, dashed lines correspond to the scattering on 
the random potential and vertexes correspond to the contact magnetic interaction Eq.2. The 
blocks of the diagrams shown in Fig. 2c correspond to the so called Cooperons and Diffusons 
P^'^i'^:^): which are tensors as functions of electron spin indexes ^^'^1. The diagrams shown 
in Fig. 2a have been considered in [^^-le] xhey give main contribution to the correlation 
function {dIij{Ri, R2)SIki{Ri, R2)) ~ l-^i — -R2I These diagrams correspond, however, to 
the correlation function {SIij{Ri, R2)SIki{R3, R4,)) which decays exponentially when |^i — 
-R3I1 I-R2 — ^4] 3> As a result, for example, in the case L ~ L2 ^ L3 ~ Li the contribution 
to {SlijSlki) from these diagrams is of the order of 

I'oinp-.n^f. (5) 

We assume the density of localized spins in the ferromagnets in of order of n. 

The diagrams Fig. 2b give much smaller contribution to < 5Iij{Ri, R2))dIki{Ri, R2) >, 
but they describe long range correlations {5Iij{Ri, R2)5Iij{R3, R4)) ~ R-'^^'^'^^ when \Ri - 
-R3I ~ I-R2 — Ra\ = R ^ I- As a result, it is these diagrams that determine the of the 
correlation function of the interlayer exchange energy {lijlki) at L ^ I. The qualitative 
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explanation of the origin of the correlation is as follows. The mesoscopic fluctuations of the 
exchange energy 5Iij{Ri, R2) result from the interference of random probability amplitudes 
of diffusion paths between the points Ri and R2. Among these paths there are some which 

— * — * 

visit points R3 and R4 (An example is the hne "a" in Fig. la). This leads to the mentioned 
above correlation between Iij{Ri,R2) and IkiiRs, Ra)- As a result, we have: 

{Slij6lki) = -llElTY^u I dRidRidR^dRi 

TV m 

{aiP:{R^, R2)akP:{R2, R^)a^P:{R^, R4)aiP:{R^, R^) + 
aiP^iRu 4)a/j(i?2, R,)akP^{R,, R4)aiP^{R4, Ri))- (6) 

Here to — 7r(2m + 1)T is the Matsubara frequency, m is an integer number T is the temper- 
ature and (Tj are spin operators. Integrations over R^ and R2, R4 in Eq.6 is performed 
over volumes of the first and the second ferromagnetic layers respectively. Results of calcu- 
lation of Eq.6 depend on the ratio between the lengths L, L2, Lt — \f^, Lgo — ^/DtTo and 
on boundary conditions for Cooperons and Diffusons. Here L^o, Tso are the spin-orbit relax- 
ation length and time, respectively, and D is the electron diffusion coefficient in the N layer. 
In the case of the "open" geometry of the N layer shown in Fig. la and Lt, Lgo ^ L > I; 
L, L2 ^ Li, L3, we have 

< 5hj5hi >= ^^^^^Xj^-{pfL,)%j6,i (7) 

3^77 2 [PFI) 

Here X is a factor, which is of the order of unity when L L2 < Lt and C,{x) is the Zeta 
function. In different limiting cases we have: 

{^Y if Lt > L2 > L 
^ iiL2>LT>L 

In the case L > Lt the expression for X acquires an additional exponentially small 
factor exp{—j^). In the case Lso > L the minimum of the exchange energy corresponds 
to parallel or antiparallel orientation of the layer's magnetizations {9 equals zero or tt). In 
the opposite limit Lgo <^ L we get the same formula as Eq.8 but without the factor dijdki- 



This means that the exchange interaction between the F layers is of the Dzialoshinski-Moria 
type and a minimum of the exchange energy corresponds to a sample specific angle 9{ni,n2) 
distributed randomly in the interval (0,7r). It is interesting that in the case L L2 < Lt 
Eq.7 turns out to be independent of L. Let us note that the diagrams shown in Fig.2a give 

While deriving the results presented above we neglected the sensitivity of the boundary 
conditions for Cooperons and Diffusons to the change of magnetization directions in F-layers. 
In the case of the open sample geometry Fig. la this is correct, provided -C 1. To get 

an estimate for {Slijdlki) in the opposite limit one has to substitute the factor A{Lipp) in 
Eq.8 for Ep. For example, in the case Lt > L L2 > Lgo we have 

{SkjShi) ~ EUpfI)-' ~ -. (9) 

T 

Here r is the elastic mean free path in the metal. We would hke to stress that the origin of 
both Eq.7 and Eq.9 is the exstance of the long range correlation of signs if the quantities 
Iij{Ri,R2) and /fc/(i?3,i?4) on distances much larger than /. 

In the case when the length of the ferromagnetic layers L2 is long enough one should take 
into account the random fluctuations of the direction of magnetization along the layers. A 
solution of this problem is beyond the scope of this paper. 

As it is usual in the physics of mesoscopic metals t^^'^^l, the external magnetic field 
changes the electron interference pattern and, thereby, 9{nin2) turns out to be a random 
sample-specific oscillating function of the magnetic field H. In the case of the open geometry 
of the N layer shown in Fig. la for L2 ~ L the characteristic period of these fluctuations 
is 5 Hi ~ ^ . Here $0 the flux quantum. There is also another characteristic 

magnetic field in the system which corresponds to the interaction energy between the external 
magnetic field and magnetic moments of the F layers which is of order of the exchange 
energy between the layers: 5H2 ~ — 3 r ^ - — 3 r^r r • Here a is the Bohr magneton. 
At large enough L, 5H2 » 5 Hi and therefore the exchange energy 7^ will be a random 
sample-specific oscillating function of the magnetic field in the situation where the interaction 
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energy between the magnetic field and the ferromagnetic moment is still negligible. As a 
result, the relative orientation of magnetizations of the F-layers, characterized by 0{ni, n^)-, 
will be a random function of the magnetic field. In the opposite limit, when 5H2 < SHi, 
9{ni,n2) monotonically decreases with the magnetic field. Even in this case it is possible to 
see the random oscillations of 9{ni,n2) near the magnetic field corresponding to spin-fiop 
transitions. 

Another way to change the relative orientations of the F-layers is demonstrated in Fig. lb. 
Namely, 9{ni, 712) is a random sample specific function of the order parameter phase differ- 
ence (xi — X2) in superconductors 5*1 and 5*2 in Fig. lb. The reason for this is that some 
diffusive paths connecting points 1 and 2 in Fig. lb can visit superconductors (line "b" in 
Fig. lb) and the corresponding amplitude of probability to travel along these paths acquire 
the additional phase {xi ~ X2) Another consequency of the phase dependence of the 
exchange energy is that the critical Josephson current of the device shown in Fig. lb de- 
pends on the angle 9 between magnetizations of F layers. At last, a change of electrical 
or chemical potentials in the metal (or in the degenerate semiconductor) also lead to the 
oscillations of lij. In the case L ~ L2 characteristic period of the oscillations as a function of 
the chemical potential /x is of the order ol5jjL ^ Ec= f^^' In the case when the degenerate 
semiconductor is a part of a MOSFET one can change fi by changing the voltage on the 
gate. 

The resistance of the system considered above is a random oscillating function of the ex- 
ternal magnetic field. There are two mechanisms which cause these oscillations. l.The usual 
mechanism for mesoscopic metallic samples: I^^'^^l amplitudes of probability to travel along 
diffusion paths acquire additional random phases, proportional to the magnetic field. The 
characteristic period of the oscillations is of the order of 5Hi. 2. The magnetic field induced 
change of 9{ni,n2) which corresponds to a change of the boundary conditions for electrons 
in nonmagnetic metal.As a result the rotation of the ferromagnetic layer's magnetizations 
leads to mesoscopic fluctuations of the resistance of the sample. The characteristic period 
of the oscillations correspond to {9{H) — 9{0)) ~ 1. 
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The calculations presented above do not take into account the dipole interaction between 
magnetic moments in different F-layers. This interaction decays with the distance between 
the layers L and can be much less than the considered above mesoscopic part of the exchange 
energy, which is independent of L if L <^ Lt- 

In conclusion, we would like to mention that the considered above mechanism can deter- 
mine random anisotropy of the exchange interaction in bulk disordered ferromagnets. Let 
us introduce an exchange field 



- f He.dT. (10) 

V Jv 



averaged over a volume v — Lq, I < Lgo < Lq < Lt- Here H^x is the local exchange 
field acting on localized spins in ferromagnet. The correlation function of this field can be 
calculated with the help of the diagram Fig. 2b. 

^HM.A?))-^.j^J(\r-r'\) (11) 

Here /(r) is the function which is equal to unity at r « Lo, decays as /(r) - {^f at 
Lq < r < Lt and exponentially decays when r > Lt- The above-mentioned long range 
correlation between SIij{Ri, R2) and 5Iki{R3, R4) manifests itself in the factor {ppLo)"^ in 
Eq.ll. An assumption about short range correlations would lead to an expression which is 
proportional to {pfLo)~^. 
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FIGURES 

Fig.l. Schematic pictures of the ferromagnet-nonmagnet layered systems. 
Fig. 2 Diagrams for the correlation function (SlijSIki). 
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Fig.l 




c) 
Fig.2 



